Orbital effect of Cooper pairs on Kondo effect in unconventional superconductors 
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A new type of Kondo effect peculiar to unconventional superconductors is studied theoretically by 
using the Wilson's numerical renormalization group method. In this case, an angular momentum 
of a Cooper pair plays an important role in the Kondo effect. It produces multichannel exchange 
couplings with a local spin. Here we focus on a p x + ip^-wave state which is a full gap system. The 
calculated impurity susceptibility shows that the local spin is almost quenched by the Kondo effect 
in the strong coupling region (7k /A — * oo), while the ground state is always a spin doublet over 
all the 7k /A region. Here 7k and A are the Kondo temperature and the superconducting energy 
gap, respectively. This is different from the s-wave pairing case where the Kondo singlet is realized 
for large 7k/ A values. The strong coupling analysis shows that the p x + ip y -w&ve Cooper pair is 
connected to the Kondo singlet via the orbital dynamics of the paired electrons, generating the spin 
of the ground state. This type of Kondo effect reflects the symmetry of the conduction electron 
system. 



I. INTRODUCTION 

The search for variety in the Kondo effect was stim- 
ulated by the discovery of the non-Fermi liquid behav- 
ior in heavy fermion materials and by the observation 
of the two-level system character in metallic glasses and 
nanoscale devices. In these cases, orbitaLdegrees of free- 
dom are very important as well as spin.EJn Due to incom- 
plete compensation of a local moment, it displays sin- 
gular behavior in physical quantities, differing from the 
ordinary single-channel Kondo effect. This Kondo prob- 
lem has been investigated mainly for metallic systems 
with no gap in the density of states. 

For a gapped system, the problem of a magnetic impu- 
rity in a s-wave-superconductor has been investigated for 
a long time.uil3 Since the s-wave property is described 
by a single channel, the Kondo effect itself is understood 
in the framework of the conventional theory based on 
the single-channel Kondo model. On the other hand, no 
one has paid attention to multichannel properties of p- 
wave and d-wave superconductors derived from angular 
momenta of their Cooper pairs. They can exhibit a new 
type of Kondo effect peculiar to the unconventional su- 
perconductivity. 

The unconventional superconductivity is very sensitive 
to impurities and surface boundaries. It displays peculiar 
phenomena not observed in the s-wave superconductiv- 
ity. It is known that scattering objects destroy the uncon- 
ventional superconductivity and induce bound statesxil 
The nuclear magnetic resonance (NMR) T^ 1 measure- 
ment performed in Zn-doped YRapCuaOv-^ showed a 
constant T\T at low temperaturesEj reflecting finite den- 
sity of states at lower energies. As shown by the low tem- 
perature scanning tunnel microscope (STM), a single Zn 
atom actually generates the low energy bound state.EJ 
For the d^^-wave, the (lr-1,0) surface is regarded as 
the pair-breaking boundaryO A very recent experiment 
showed a zero energy peak in the tunneling conductance 



between the (AJ.,0) boundary of YBa2Cu30y_5 and a 
normal metal.E^I Thus much attention has been paid to 
the effects of non-magnetic impurities and static bound- 
aries in the last decade. 

Recently the effect of magnetic impurities in uncon- 
ventional superconductors has attracted attention grad- 
ually with the development of experimental studies. Ac- 
cording to the NMR experiments, Il3 a local moment is 
induced around a Zn site substituted for a Cu site in 
YBa2Cu307_,5. It is expected that, the observed Kondo 
behavior is due to this moment .113 Similar behavior is 
also observed in Li-substituted YBa2Cu307_5.E3 In this 
case, the local magnetic susceptibility of the Li induced 
moment shows strong reduction of the Kondo screening 
below the superconducting transition temperature. At 
present the search for the influence of unconventional su- 
perconductivity on the Kondo effect is still in progress. 
We are also interested in the Kondo effect for differ- 
ent types of superconductors other than such high-T c 
cuprates. 

Several studies of a magnetic impurity in a conven- 
tional (s-wave) superconductor clarified how a bound 
state generated around the impurity depends on the 
Kondo_temperature Tk and the superconducting energy 
gap A.ETtJ The numerical renormalization group (NRG) 
study by Satori et aln showed that the bound state en- 
ergy is scaled by Tk/ A and the interchange of singlet and 
doublet ground states occurs smoothly. Even if an energy 
gap exists in the density of quasiparticle states, a local 
spin couples with the quasiparticlcs antiferromagnetically 
to form a spin singlet (Kondo singlet) for Ik /A > 0.3. A 
similar competition was discussed in gapless.(i :z .2_j / 2-wave 
superconductors such as high T c cuprates.lia In this case, 
the gapless density of states stabilizes the Kondo singlet 
more favorably than in the gapped s-wave case. However, 
it is not certain that a singlet ground state is always re- 
alized in all types of superconductors when the Kondo 
effect is strong enough. A Cooper pair in the unconven- 
tional superconductivity has an angular momentum and 
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can produce multichannel scattering of quasiparticles at 
an impurity. As well as for metals, a difference between 
single channel and multichannel properties in the Kondo 
effect can be expected for superconductors. 

In order to elucidate the multichannel property, we 
study the Kondo effect due to a local spin coupled 
to a p x + ip y (or d x 2_ y s + id xy )-wave superconducting 
state. The c£-vector for the p x + ip y -~w&ve is given by 
d k = z(k x + iky) where z is a unit vector in the z axis. 
Since the density of states has a full gap like the s-wave 
case, we can find the peculiarity of the unconventional 
superconductivity in comparison with the s-wave. The 
property of the p x + ip v superconductivity is analogous 
to that of the A phase found in the superfluid 3 He. A 
lot of attention has been paid to this superconductiv- 
ity since the .discovery of the spin triplet superconductor 
Sr 2 Ru04.IErc3 The recent theoretical studies proposed a 
number of unique physics associated with the brokenlime 
reversal symmetry observed by the /iSR experiment .E3 In 
the p x + ipy-w&ve superconductor, spontaneous current 
appears close to the s catte ring objects and it induces 
a finite magnetic field.oS This argument explains the 
/iSR data wellBa If the system induces the spontaneous 
magnetic [field by itself, a spontaneous Hall effect can be 
expected£3 A vortex charging effect is. also interesting, 
peculiar to the rotating Cooper pairsH3 

The most important point in this paper is that for the 
Px + ipy-w&ve case one of the Cooper pairing electrons is 
coupled to a local spin and the other has no direct cou- 
pling with the local spin. As we will discuss later, the lat- 
ter can interact with the local spin through the supefeptu 
ducting order parameter. We use the NRG methodEHTLil 
to study this new type of two-channel Kondo effect. As 
for the s-wave case, we discretize the bare conduction 
band logarithmically and apply the Bogoliubov transfor- 
mation to obtain a diagonal form of NRG Hamiltonian 
for the superconductivity. In the NRG Hamiltonian, the 
two independent channels consisting of the I = and 
/ = 1 orbitals are coupled to each other at the impurity. 
The NRG result shows that the ground state is a spin 
doublet over all the Tk/A region and Ik/ A dependence 
of a bound state expresses a competition between the 
Kondo effect and the superconducting energy gap. The 
solution of the strong coupling limit (Tk/A — > oo) shows 
that the ground state is a spin doublet due to dynamical 
coupling between the I = and I = 1 electrons, while the 
local spin is almost quenched by the I = electrons. 

This paper is organized as follows. In Sec. II we de- 
rive the NRG Hamiltonian to investigate the two-channel 
Kondo effect specific to the p x + ip y (or d x 2^ y 2 + id xy )- H = B (i?kin + H A + H lr 
wave superconductivity. Here the origin of the two chan- 
nels is explained thoroughly. In Sec. Ill we show the 
NRG results in comparison with the s-wave case. We 
also calculate a local spin susceptibility and discuss the 
Ising case. It is very useful to understand the quantum 
effect of the magnetic impurity in the gapped electron 
system. This paper is closed with concluding remarks in 
Sec. IV. 



II. FORMULATION 

In this section we derive the NRG Hamiltonian for the 
Kondo effect in the p x + ip y -w&ve (d x 2_ y 2 + id xy -wave) 
superconductor. Here we study a single magnetic im- 
purity at the center in two-dimensional superconducting 
systems. 



A. Model Hamiltonian 

Let us begin with the following Hamiltonian: 
H = H kin + H A + H imp , (2.1a) 

-ffkin = ^2 £ fcaL a ^' ( 2 ' lb ) 
krr 

H ^ = E ( A fc a fc T a -fci + H - c -) > ( 2 - lc ) 

fc 

#im P = - E ■ cr aa , +V6 a ^)al, a a k ' a ,. (2. Id) 

kk'crcr' 

Here k and a represent the wave vector and spin of 
the conduction electrons, respectively. a ka is the anni- 
hilation operator of the conduction electron, e k is the 
kinetic energy of the conduction electron, and Afc is the 
momentum dependent order parameter. S expresses the 
5 = 1/2 spin operator for the magnetic impurity, and <x 
is the Pauli matrix for the conduction electrons. J(< 0) 
and V are the antiferromagnetic and non-magnetic cou- 
plings, respectively. Here we have assumed a short-range 
(s-wave) scattering impurity. For the p x + ip y -w&ve pair- 
ing state, the order parameter has a form of A k = Ae'^fc . 
Here A is the order parameter which we assume a real 
number, and 4>k is the angle of the wave vector measured 
from the k x axis. Since the orbital angular momentum 
is a good quantum number for the p x + ip y -wave, it is 
convenient to transform the wave vector representation 
into the following two-dimensional polar coordinate rep- 
resentation: 



a ka = Y,(^) l \Ji+i{kR)\^a klal 



(2.2) 



where Ji is the l-th Bessel function, and R is the radius of 
the two-dimensional system. I is the z-component of the 
orbital angular momentum of the conduction electron. 



Substituting Eq. (2.2) into Eq. (2.1), we obtain 



#kin = j dkk^2 a\ lcr a kla , 

^~ la 

B A = f dkYj n ; ; (iAo£ lt «,,_, +u 
I 

Hi m p = (-JS ■ <J aa i + V8 a ^:) f 00cr foOa' 



(2.3a) 
(2.3b) 

(2.3c) 
(2.3d) 
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— 7= J dka k ,i= ,, 



(2.3e) 



_where B = k F R/(nm) = 2k F RE F /Tr, A = A/(2E F ), 
J = Jp, V = Vp. Here E F = k F / (2m) is the Fermi 
energy, and p = mR 2 /2 = (k F R) 2 / (AE F ) is the density 
of states at the Fermi energy for the two-dimensional 
conduction band. In Eq. J2.3| ) we have used the follow- 
ing asymptotic form of the Bessel function: Ji(kR) ~ 
J2/(nkR) (for a large kR). The wave number k in Eq. 
(2.3) is normalized by k F . The kinetic energy of the con- 
duction electron is linearized near the Fermi energy. For 
simplicity, we have used the same values for both the su- 
perconducting cutoff energy and the band width fcpinSince 
it is assured that this does not alter the resultsfl The 
fermion represented by the operator /oo CT corresponds to 
a fermion localized at the impurity. The first subscript 
of /ooo- represents the 0-th shell of the NRG Hamiltonian 
which will be derived later (see Appendix). The second 
zero means the I = orbital. Notice that the angular mo- 
mentum of the p x + ipy-w&ve Cooper pair equals to one 
(see Ha), and that the I = and I = 1 orbitals are cou- 
pled. Since the angular momentum is a good quantum 
number for the p x +ip y -w&ve, both orbitals are decoupled 
from the other orbitals. Due to the short-range (s-wave) 
scattering, the electrons only in the / = orbital couple 
with the impurity. The electrons in the I — 1 orbital can 
also interact with the impurity via the p x + ip y -wave or- 
der parameter. It is suffici ent to choose only the I = 
and I — 1 orbitals in Eq. (2.3) to study the Kondo ef- 
fect for the p x + ipy-w&ve state. Therefore we can apply 
the NRG method to this Kon d o p roblem as we do to the 
two-channel Kondo problem. HEj We rewrite then -ffkin 
and Ha as 

Hkin = J dkk^ { a k0a a k0a + a ll<x a ^) ' ( 2 ' 4& ) 

H A = J^dkJ2(iHiA ,-*+K-' 



(2.4b) 



We can obtain a similar Ha for a d x 2_ y 2 + id xy -wave if 
we replace iA by — a A and / 
order parameter for the d x i_ %j i 
A fc = Ae l2<#, fc. 



= 1 by I = 2, since the 
■ id xy -w&ve has a form of 



B. Transformation of the Hamiltonian 

In .our formulation we follow the procedure by Sakai 
et alB and derive a NRG Hamiltonian. There is another 
way to derive the Hamiltonianja which we describe in Ap- 
pendix. First we diagonalize the conduction electron part 
by introducing the following Bogoliubov transformation: 



a kla =u k u k _^+v* k a\ + 



(2.5a) 
(2.5b) 



Here a k T -± a is the annihilation operator for the super- 
conducting quasiparticle with a spin in r channel. The 
coefficients u k and v k are defined by 




/ 1 \k\ 



(2.6a) 
(2.6b) 



where E k = (k 2 + A 2 ) 1 / 2 is the energy of the quasipar- 
ticles. The p x + ip y -w&ve Hamiltonian is then expressed 
in terms of the superconducting quasiparticle as 



H, 



-ffkin + Ha 



I dkE k Sgll(k) a kra a kra- (2-7) 

For the later convenience, we introduce the following 
transformations: 



^|fc|,-T,- 



'\k\ 



U-\k\ 



~j= ( a kra ~ iT fikro) : 



B kra = ( a kra + ir Pkra) , 

C+, CT (fc) = A ki+<a , 

C+A-k) = -iBl,-,-*> 
C-,cr{k) = iA k __ a , 
C-, a {-k) = B k .„. 



(2.8a 
(2.8b 

(2.8c 

(2.8d 
(2.8e 
(2.8f 
(2.8g 



The Hamiltonian ( p.7j ) can be expressed by the above C 
operators as 

H Px+iPy = f dkE k V TsCl a (sk)C Ta {sk). (2.9) 

JO „ L „ „ 



For the impurity part, the operator /ooo- is also ex- 
pressed in terms of the superconducting quasiparticle as 



dk 



V2Jo 



u k (a k+a + $l__ a ) 

w k+ (A k+a + iAl__ a ) 
w k _(B k+a -iBl__ a ) 



(2.10a) 
(2.10b) 



By introducing the following operators 
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a T ,a = N, 1 dk w k+ A k 



b T ,a = N_ x dk w k ^B ktTi(7 , 
Jo 



N± = / dk w 



1/2 



(2.11a) 
(2.11b) 

(2.11c) 



foou is expressed as 



foOa 



1 f 



V2 



N+a, 



iN-V 



+ i(N+a^_ a -iN-b^ cr ) 



1 

71 



(C0, + ,(J + C0,_, CT ) , 



(2.12) 



where co. r ,o- is defined by 



C0,-,a = I 



J dk (w k+ A k + a - iw k _B ] k _ 
/ dk V] u> fc . s C+, CT (sA;), 



(2.13a) 



= y dk [w k+ iAl _ _ a + w k _B k + i 
= [ dky^w ks C-. a {sk). (2.13b) 

J s=± 

The exchange Hamiltonian Hi mp is then written as 

-ffimp = - (-JS- CT aa ' + VS aiCr /) cJ TCT C 0r / o -/ . 

(2.14) 



In Eq. (2.14), the subscript r = ± represents the two 
channels, both of which consists of the I = and I = 1 
orbitals. While t he im purity interacts with only the I = 
electron, Eq. ( 2.14 ) contains the channel-flip (X} TT ') 
term s. Th is is because the operator /ooo- is e xpress ed as 
Eq. (2.12). The two-channel property in Eq. (^T4|) orig- 
inates from the angular momentum of the p x + ipy-w&ve 
Cooper pair. In order to treat the Kondo problem in 
the superconducting state, first we have diagonalized the 
Hamiltonian of the conduction electron part. The diag- 
onalized representation (r representation) in this Hamil- 
tonian is given by the linear combination of the I = 
and I — 1 orbitals, since the p x + zp^-wave mixes the 
two orbitals. It is natural that the impurity part of the 
Hamiltonian contains the channel-flip part in the r rep- 
resentation. 

In contrast to the present p x + ip y -wetve case, the stan- 
dard two-channel Kondo model was derived frana-the in- 
ternal degrees of freedom of the impurity atom.ErE-3 In our 



model, the impurity does not have such internal degrees 
of freedom except for the S — 1/2 spin. However, we can 
obtain the two-channel property due to the orbital effect 
of the Cooper pairs. Therefore the present two-channel 
property for the p x + ip y -w&ve can provide new Kondo 
physics peculiar to the unconventional superconductors. 



C. Discretization 



Following Wilson's procedure)^ we discretize the con- 
duction band with a logarithmic mesh in the k space. 
The Hamiltonian for the discretized conduction band is 
written as 

oo 

Hp*+iPy = Y ^2 TsE mCl a (sm)C Ta (sm), (2.15a) 



m— res 



E m = ^kl + A^, 

fc m = (l+A- 1 )A- m /2, 



(2.15b) 
(2.15c) 



where A > 1 is the discretization parameter. The above 
C operators satisfies the Fermion commutation relation: 

|(7 T(T (sm),C'|, cr ,(s'm')| = 8 T . T >5 a ^5 s , s >5 m . m > . (2.16) 

The fermion operator for the localized orbit (the 0-th 
shell orbit) cq t<7 is also expressed with the C operators 
as 



Co 



ra = (1 - A" 1 )" 1 / 2 Y, w ms K- m l 2 C Ta (sm), (2.17a) 



W ms = J (1 + S ^-)/ 2 - 



(2.17b) 



We define c nr(T for the n-th shell orbit and satisfies the 
following commutation relation: 

\pnra 7 Hp x +ip y ] — ^n— l^n— l.rcr "1" TTj n C nT(7 

+ tnCn+l.TO- (2-18) 

For n = we have 

[co T!7 ,H Px+iPy \ =(i-A-y/ 2 

x Y w ms A~ m/2 TsE m C T(T (sm) 

sm 

= Tr/oCOrcr + t C lTa , (2-19) 

where 770, io, and c\ ra are defined by 

Vo = (l-A- 1 )J2™ 2 ms^~ m sE m , (2.20a) 



to 



(l-A-^^LA-"^™-^) 2 

sm 

(l_ A -l)l/2 



1/2 



(2.20b) 



to 
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w sm A m/2 (sE m - T] )C Ta (sm). 



(2.20c) 



Repeating this procedure, we obtain the total Hamilto- 
nian which-has the following hopping type with staggered 
potentials n 



(2.21a) 



H 



(2.21b) 



Hi mp = - ^2 {-JS-cr arT i +V6 a<a ') 4 to .c 0t v- (2. 



21c) 



Here ry n and t n are given by 
r, n = (-1)"A, 

e« = [l - A-(" +1 ) 
x "i_A-( 2 " +3 > 



1 - A~ (2,l+1) 

-1/2 



-1/2 



(2.22a) 
(2.22b) 

(2.22c) 



In the practical calculation, we use the following scaled 
Hamiltonian in a recursion relation given by 



H N+1 = A 1/2 H N + ^ [eN(4f +ltT(T c NTa + H.c.) 

rcr 

+ (-l) N A N / 2 rAJ N+hTa c N+hTa 

rr'crcr' 

-^tAc^co^Ja- 1 / 2 . 

rcr 

Here the parameters are defined by 
A= A 



(2.23a) 



(2.23b) 



J = 
V = 



(l + A-i)/2' 
J 

(l + A-i)/2' 
V 



The Hamiltonian ( |2.2l| ) is related to Hn and is given by 
1 + A" 1 



H 



lim A-^- 1 )/ 2 ^ 



JV- 



(2.25) 



In Eq. ( 2.2S ) , cjsstc is an operator of the NRG fermion 
quasiparticle in the iV-th shell. The subscript r = ± 
represents the two channels which are constructed by the 
/ = and / = 1 orbitals as discussed in the last part of 
the previous subsection. We can obtain the same NRG 
Hamiltonian for a d x _ 2 -y 2 +id xy -w&ve. In order to check 
reliability of numerical results, we compare different data 



by changing cutoff of the eigenstates and the discretiza- 
tion parameter A. In the NRG results we show through- 
out this paper, we keep the lowest-lying ~ 500 states at 
each renormalization step and take A — 3. For calculat- 
ing a local spin susceptibility, we keep <~ 1200 states to 
reduce the cutoff dependence. 



III. KONDO EFFECT IN P-WAVE 
SUPERCONDUCTOR 

In this section we present our results for the Kondo 
effect in the p x + ip y -w&ve superconducting state. We 
show numerical results and then discuss the strong cou- 
pling limit to see how the orbital dynamics of the Cooper 
pairs affects the Kondo singlet. 



A. Numerical results 

First we discuss a A = case where the only I = 
electrons couple with the local spin and the I = 1 elec- 
trons are free from the impurity. The NRG result shows 
that the Kondo singlet is realized by the relevant I = 
electrons at the stable fixed point. This fixed point is of 
strong coupling type (\J\ — > 00) and the fixed-point NRG 
energy levels is reproduced by truncating the N = site 
(impurity_aite) from the I — chain of the conduction 
electronsrj On the other hand, the fixed point for I = 1 
is of free-electron type (J = 0). One of the two chains 
produces a zero energy state in one-particle energy levels 
at the fixed point. We can fill at most two particles (spin 
up and down) into the zero energy state and we have 
quartet ground state for A = 0. 

Next we discuss a J = V — case with finite A. In 



this case, the NRG Hamiltonian (2.23) consists of the two 
independent r = ± channels related to the superconduct- 
ing state. The Hamiltonian of each channel is same as 
that of the free s-wave superconductivityJS At the odd 
renormalization steps, one-particle energy level for the 
free s-wave is particle-hole symmetric. At the even steps, 
there is an additional one-particle energy level which is 
hole-like (particle- like) for a positive (negative) A. In 
the present p x + ip y -w&ve case, one particle energy level 
is symmetric at all the renormalization steps, since the 
signs of A for the r = ± channels are different each other. 
In one-particle energy levels, there is no zero energy state 
due to the finite A. The ground state is obtained by fill- 
ing particles into the energy levels and we have a single 
ground state. 

Let us turn to the Kondo problem in the p x + ip y -~w&ve 
state. In this case both A and J are finite. Since J grows 
up as k N / 2 with increasing N, it is relevant in the Kondo 
effect . In a ddition, A is also relevant since it grows as in 
Eq. ( [2.23 ). Therefore we can expect a competition, be- 
tween the Kondo effect and the superconductivity!! As 
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FIG. 1. Tk/A dependence of the bound state energy 
(Eb/A)* for p x + ipy-wave. The ground state is always a 
spin doublet and the bound state level is measured from the 
ground state for each Tk/A. The first excited state is a par- 
ticle-hole doublet with no spin for V = 0. The circle and star 
represent the result with fixed Tk for Tk = 9.20 x 10~ 5 and 
1.76 x 10 -3 , respectively. The triangle-up and triangle-down 
are the results for fixed A = 0.001 and 0.005, respectively. 
Here A = 1.5A. For both Tk and A, a unit of the energy is 
the band width 2Ep. 



the renormalization step increases, the excited state en- 
ergies increase with An ~ A N / 2 , which is the energy gap 
at the N-th renormalization step. Since the magnetic im- 
purity can destroy the superconductivity, we have bound 
states below A at. The ratio of the bound state energy 
and An approaches a constant value with the increase 
of the renormalization step. We express the convergent 
value of the ratio by (E B /A)*. Figure shows the T K /A 
dependence of the first excited state (bound state) for 
V = 0, where the Kondo temperature is denned by 



-i/VD- 



(3.1) 



Here, a unit of the energy is taken to be the band width 
2Ep ■ Wc find that the bound state energy level is scaled 
by Tk/A as in the s-wave case. In a small Tk/A region, 
the ground state is a spin doublet as expected, and the 
first excited state is a particle-hole doublet with no spin. 
With the increase of Tk/A, the energy level of the bound 
state (the first excited state) decreases (see Fig. [j]). As 
Tk/A — > oo, it approaches the ground state energy and 
the ground state is four-fold degenerate. The result for 
Tk/A — * oo is connected smoothly to that for A = 0. 

The most important point in Fig. [j] is that the ground 
state is always a spin doublet. This means that the spin 
of the ground state survives even in the strong coupling 
limit for the p x + ip y -wa.ve superconductivity. In the s- 
wave case, the Kondo singlet is realized in the large Tk/A 
regionO, resulting in a spin singlet ground state. Thus the 
Kondo effect for the p x + ip y -wave is quite different from 
that for the conventional s-wave. 

While the ground state is always a spin doublet for 
the p x + ipy-w&ve, we find that the Tk/A dependence 
of the bound state energy level is very similar to that of 
the s-wave. For the s-wave, the Kondo singlet is real- 



ized in Tk/A > 0.3 region against the superconducting 
energy gapE In Fig. (jl the curvature of the (Eb/A)* 
curve found around (Eb/A)* = 0.5 is very similar to 
that of the s-wave case found with the interchange of 
the ground state. This represents a crossover from the 
A dominant region to the Tk dominant region. It oc- 
curs around Tk/A = 0.3 as in the s-wave case. Thus 
the Kondo effect can overcome the energy gap also in the 
Px + ipy-w&ve case, although the ground state is still a 
spin doublet. This means that the competition between 
the Kondo effect and the energy gap is characterized by 
such Tk/A dependence. All the above results hold for 
the d x 2_ y 2 + id xy -wave as well, since we can obtain the 
same NRG Hamiltonian. 

Finally, we me ntion the effect of the potential scat- 
tering V in Eq. (2.23). It is well known that a single 
non-magnetic impurity destroys the unconventional su- 
perconductivity. It can generate a bound state locally 
around the impurity. If we treat the magnetic impu- 
rity realistically, we have to start from the Anderson 
model. Since the generalized Kondo model consists of 
both magnetic and non-magnetic interactions, we have 
to take both couplings into account. Since the effect of 
V breaks the particle-hole symmetry, it lifts the degen- 
eracy of the first excited particle-hole doublet. However, 
the ground state is still a spin doublet and V is irrelevant 
even if it is large enough. 



B. Strong coupling limit 

Here we study the mechanism of the spin doublet 
ground state for the p x + ip y -wave. Let us discuss the 
strong coup ling l imit (\J\ — > oo) case. The Hamiltonian 
Hq in Eq. (2.23) can exhibit the solution in this limit. 
Returning the channel (r) representation to the angular 
momentum (I) one, we can rewrite the Hq for V = as 



Ho = Hq j + H a , 
H j = -JJ2 S - (t 



crcr' 90cr90cr' , 



H 0A = -A^2(gl a g 1(T + g\ a goa), 

where the following transformation has been used: 
1 



cn=o.- 



V2 



(gi=o,a + Tgi=i l<r ) . 



(3.2) 



(3.3) 



Here gi=o,a and gi=i, a correspond to / 00ct and /J Xo ., re- 
spectively. Notice that the order parameter corresponds 
to the hopping parameter between the I = and / = 1 
local orbital sites, since the p x + ip y -w&ve Cooper pair 
consists of the I = and I — 1 particles. When \J\ — > oo, 
the I = particle couples with the local spin strongly to 
form a spin singlet |s): 



|B> = ^(flStU>-flSiit>). 



(3.4) 



G 



Ipl'} 




lpl> 

(a) (b) 
FIG. 2. Excited states generated from the unperturbed 
ground ones by the perturbation of A. (a) Particle-hole dou- 
blet case, (b) Spin doublet case. 



local spin 



l Cooper pair 




I = 



I = 1 



(a) 




Kondo 
singlet 



I = 



I = 1 



(b) 



where | f) and | J.) represent the local spin states. Since 
the hopping is forbidden, the ground states are 4-fold 
degenerate for \J\ — > oo: 



|s T> = ffi T |s), 




l«l)=.9 li |s), 








|p2) = s} T s}js), 


(3.5) 



where |scr) (a =|, |) and |pi) (i = 1, 2) are the spin dou- 
blet and particl e-ho le doublet, respectively. At a finite 
| J |, A in Eq. (3.2) lifts the degeneracy. The pertur- 
bation (the A term) ge ner ates excited states from the 
ground states in Eq. (|3.5| ). We show the connection 
among these states in Fig. 0. One of the particle-hole 
doublet states is connected to a single excited state. For 
|pl), the excited state is given by [see Fig. ^ (a)]: 



1 



|pi') = ^(^ T |4)-^JT) 



(3.6) 



The similar argument holds for |p2). We can obtain the 
following eigenvalue for the particle-hole doublet up to 
the fourth order of the perturbation: 



3 2 / A \ 8 / A 

2 H^V + ^{W\ 



\J\. (3.7) 



For one of the spin doublet states, it is connected to 
four excited states by the perturbation up to the fourth 
order. For |s f), these excited states are given as follows 
[see Fig. |(b)]: 

|0T) = <MJT), 
I 1 T) =£1^11 1 T). 

|tl) = .9u^ T l T>, 



|t2) = 7f^ T [AU)+foilt) 



(3.8) 



and the similar connection is obtained for |s J.) as well. 
The energy eigenvalue can be obtained from the following 
equation: 



FIG. 3. Schematic picture of the ground state wavefunc- 
tion. (a) Weak coupling limit (7k /A — > 0). (b) Strong cou- 
pling limit (T K /A -» oo). 

e(e - 1/2) [e 3 + e 2 - (Ax 2 + 3/4)e - 4a; 2 ] = 0, 

e = E/\J\, x=\A\/\J\. " (3.9) 

After expanding the lowest e in x up to the fourth order, 
we obtain the energy for the spin doublet as 



Es=l/2 





\J\. (3.10) 



The energy difference between the two doublets is given 
by 



2 - A 

Es=o - E s= i/ 2 = -\J\ I j 



(3.11) 



Thus, the order parameter A stabilizes the spin doublet 
even if it is very small. 

In the s-wave case, the superconducting part of Hq is 
given by 



H, 



OA 



-A 



1 



(3.12) 



We notice that it is not given by hopping between the 
sites but is given by a potential on the / = site. This 
is because the s-wave Cooper pair is formed by only the 
I = quasiparticles. In the |J| — > oo limit, we have a 
Kondo singlet ground state as expected. 

Comparing the p x + zp^-wave result with the s-wave 
one, we conclude that the spin of the ground state is 
realized by the orbital dynamics of the p x + ip y -~w&ve 
Cooper pairs. This point is illustrated by Fig. ||. In 
a weak coupling region, the spin of the ground state is 
mainly produced by the local spin itself [see Fig. |](a)]. 
In the strong coupling limit, however, a Kondo singlet 
is formed by the I — conduction electrons. Therefore 
the local spin is almost quenched by the I = electrons. 
In this limit, one of the paired electrons having / = 1 
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FIG. 4. (a) Temperature dependence of T^iocai for vari- 
ous A. Here A = 1.5A. The exchange coupling is fixed as 
J = -0.25 (T K = 5.34 x 10" 4 ). The unit is (so^b) 2 , where 
go is a bare p-factor of the local spin, (b) Tk/A dependence 
of the effective g-factor. The circle, star, and triangle-up rep- 
resent the results with fixed T K = 9.20 x 10~ 5 , 5.34 x 10~ 4 , 
and 1.76 x 10 -3 , respectively. For both Tk and A, a unit of 
the energy is the band width 2_Ep. 



angular momentum is connected weakly with the Kondo 
singlet to gain the superconducting condensation energy. 
The latter electrons generate the ground state spin as 
shown in Fig. ||(b). 




FIG. 5. Bound state energy calculated by the NRG 
method for the Ising case. It is normalized by A. The solid 
line is the analytic solution given by Eb/A = [l + (aj z ) 2 ]~ 1,/2 , 
where J z = 1.5J Z , and a is a constant. We have used 
a = 0.465 to fit the NRG result. 



H ^ Ha-gnnSzhk-V 2 . 



(3.15) 



Therefore xi OC ai represents magnetic response of the local 
spin. In Fig. |](a) we show the temperature dependence 
of xiocai for various A with a fixed exchange coupling 
J. We can see that Txi oca i decreases monotonically as 
we decrease temperature (increase the renormalization 
step) for each A, and that it becomes constant at low 
temperatures, indicating the Curie law. However, the 
effective g-factor of the local spin depends on both J and 
A shown in Fig. |](b). We can see that g is scaled by 
Tk/A as in the case of the energy of the first excited 
state [see Fig. This result indicates that the local 
spin is actually screened by the Kondo effect, although 
the ground state is a spin doublet. 

In order to see the quantum effect (Kondo effect) of 
the local spin, we compare the result of the bound state 
energy for a Ising (classical) spin [see Fig. 0). The bound 
state energy position is given by 



C. Impurity susceptibility 

Next we see how the spin of the magnetic impurity 
is screened by the Kondo effect. For this purpose we 
calculate a local spin susceptibility which is defined by 



Xiocai(r) = gfi B 



(S z 



\h->0) 



(S z 



Tr [S z eM~PH N )] 



Here S z is the local spin operator, and (3 ~ 1 is, 
the calculation. The temperature T is given bj 

T = -(l + A- 1 )A-( JV - 1 )/ 2 / J 5. 



(3.13a) 
(3.13b) 



(3.14) 



A small magnetic field h is applied only at the local spin 
s ite. The impurity part of the Hamiltonian Ho in Eq. 
(2.23b) is replaced as 



En = 



ttJVo, 



4ft 



\Jz 



(3.16a) 
(3.16b) 



which is e^entially same as in the non-magnetic impu- 



rity case.KEI Here No is the density of state at the Fermi 
energy in the normal state, fl is the volume of the sys- 
tem. The bound state energy decreases as we increase 
the Ising coupling \J Z \. This point is similar to the result 
of the quantum spin case shown in Fig. |l|. However, the 
bound state energy is independent of A for the Ising spin, 
while it is a function of Tk/A for the quantum spin. In 
addition, we have to use an extremely larger |J 2 | value 
(\J Z \ ^S> 1) for the Ising spin to obtain a small (.Eb/A.)*. 
This indicates that the spin-flip effect (quantum effect) is 
essential to our result, and that the Kondo effect actually 
occurs in the p x + ip y -w&ve superconductor. 
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IV. CONCLUSION AND DISCUSSION 

We have studied the novel Kondo effect due to the 
orbital effect of the Cooper pairs of unconventional su- 
perconductors. In particular, we have focused on the 
Px + iPy-wave (d x i_ y i + id xy -wa,ve) pairing state. Regard- 
ing the magnetic impurity, we have assumed an s-wave 
(short-range type) impurity scattering, which is the sim- 
plest but is primarily important in real materials. The 
characteristic point of the unconventional superconduc- 
tor is that the Cooer pair has orbital degrees of freedom. 
Therefore the superconducting quasiparticles consist of 
electrons with various angular momenta. The object of 
this paper is to clarify how the orbital dynamics of the 
Cooper pairs affects the Kondo effect in unconventional 
superconductors. The followings are our main results. 

(1) p x + ipy-w&ve 

In the superconducting state, the low-energy excita- 
tions are moved out of the superconducting energy gap. 
The p x + ipy-w&ve superconductor is a fully gapped sys- 
tem. Unlike the s-wave superconductor, the p x +ip y -w&ve 
Cooper pairs have the orbital degrees of freedom. The lo- 
cal spin couples with only I = electrons which are paired 
with the I — 1 ones. However, the I — 1 electrons can in- 
teract with the local spin via the superconducting order 
parameter. We note that there is no direct coupling be- 
tween the local spin and the I = 1 electrons. In contrast 
to the two-channel Kondo. effect derived from the atomic 
structure of the impurityJII the present two-channel prop- 
erty is the consequence of the internal degrees of free- 
dom of the conduction electron system. Our NRG result 
shows that the ground state is a spin doublet over all 
the Tk/A region. This is completely different from the 
result of the s-wave superconductor in which the Kondo 
singlet is realized in a large Tk/A region. We also find 
that the first excited energies are scaled by Tk/A. In the 
strong coupling limit (Tk/A — > oo), it is found that the 
effective g-factor of the local spin is strongly reduced by 
the Kondo effect, while the p x + ip y -w&ve Cooper pair 
couples with the Kondo singlet, generating the spin of 
the ground state. Thus the orbital effect of the Cooper 
pair is important in the Kondo effect we have discussed 
here. 

(2) Origins of the multichannel property 

Generally, there are two origins of the multichannel 
property in the Kondo effect. The first one is the internal 
degrees of freedom of the impurity atom as pointed out by 
Nozieres and Blandin.EJ In this case, the channels of the 
conduction electron system is controlled by the atomic 
structure of the impurity atom. 

The second one is related to the conduction electron 
systems. In this paper we have investigated the Kondo ef- 
fect in the unconventional superconductors. As shown in 
the p x + ipy-w&ve case, the multichannel property comes 
from the angular momenta of the Cooper pairs. The in- 
dependent channels of the superconducting quasiparticles 
couple with each other at the impurity even if there is no 



internal degrees of freedom in the impurity atom (just a 
local spin). This is because the symmetry of the conduc- 
tion electron system is different from that of the impurity 
scattering. As a different variety besides the supercon- 
ducting systems, we can introduce strong momentum de- 
pendence into the bare conduction band. For example, a 
singular pseudo-gap in the vicinity._af the Fermi energy 
produces unusual Kondo behaviorO Thus the internal 
degrees of freedom of both the impurity and the conduc- 
tion electron system are important for the Kondo effect 
in real materials. 
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APPENDIX A: ANOTHER DERIVATION OF 
THE NRG HAMILTONIAN 

In Appendix, we pre sent another derivation_of the 
NRG Hamiltonian ( [2.25 ), following Satori et al.B First, 
we apply the logarithmic discretization to the conduction 
band. Then we obtain the hopping type of Hamiltonian, 



#kin = 



1 + A- 



x EE rn/3£ »(4/»+y.+/! +1 ,j^)- (ai) 

n— la 

for the free electrons. Here /.L and /„; CT are fermion 
operators corresponding to the conduction electrons with 
angular momentum I. For the BCS pairing interaction in 
the p x + ipy-w&ve superconductors, we can also obtain 
the following Hamiltonian: 

oo 

HA = E(- iA )(/n,0,<x4,l,-a ~ /n,l,-<xiW). (A2) 



To diagonalize Ha, we use the following Bogoliubov 
transformation: 



^n,+,<T — /^(/njO.CT (t)> 
bn,—,—u = —J={fn,l,-a — */n.o,o-)i 

and we obtain 

n— \ to / 



(A3a) 
(A3b) 

(A4) 



The impurity Hamiltonian is expressed by the Bogoli- 
ubov operators 6jj T(T and b nTa as 
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-ffimp — — J J]] S ■ <J aa i /oOcr/oOcr' 



ijij' 

x (&0,4 



i6A 



(A5) 



The Hamiltonian for the conduction electrons is given by 

n— tr 



x (&n, + .cr&n+l,-,-cr 



u n, + ,a u n+l,- 
+l,+,a°n, 



b n ,-,-(r + ^n+l, 



(A6) 



Next we apply the following particle-hole transformation: 

C \n, + ,a = ^2n,+,CT' 4n,-,er = *&2n,-,-trj 



-"2n-l,+,<r 



'2n-l,- 



= -6Jn-l,+,«r. (A7) 



In this manner we can exchange the roles of particles 
and holes to keep the number of fermion quasiparticles. 
Then we obtain the same form of the Hamiltonian given 
in ( p. 23 ) . The above argument can be applied to the case 



of d x 2_ y 2 + idxy-w&ve superconductors by replacing iA 
to — erA. In this case, the Bogoliubov transformation is 
given by 



4,+,<r — — /=(/l,0,<7 + a fn,2,-a), 



1 



V2 



(/n,2,-<T - f/n.Ccr)) 



(A8a) 
(A8b) 



and the roles of particles and holes are exchanged by the 
following transformation: 



C 2n, + ,a ~ ^2n,+,CTi C 2n,- .t 



-(Xbln.-.-n 



i-l,+,<7 = -^2n-l - -<r, 4n-l,-,<r = 4n-l,+,ff- ( A9 ) 
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